Electrostatics of Edge States of Quantum Hall Systems with Constrictions: Metal-Insulator 

Transition Tuned by External Gates 
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The nature of a metal-insulator transition tuned by external gates in quantum Hall (QH) systems with point 
constrictions at integer bulk filling, as reported in recent experiments [S. Rodaro et al, Phys. Rev. Lett. 97, 
046801 (2005)], is addressed. We are particularly concerned here with the insulating behavior — the phenomena 
of backscattering enhancement induced at high gate voltages. Electrostatics calculations for QH systems with 
split gates performed here show that observations are not a consequence of interedge interactions near the point 
contact. We attribute the phenomena of backscattering enhancement to a splitting of the integer edge into 
conducting and insulating stripes, which enable the occurrence of the more relevant backscattering processes of 
fractionally charged quasiparticles at the point contact. For the values of the parameters used in the experiments 
we find that the conducting channels are widely separated by the insulating stripes and that their presence 
alters significantly the low-energy dynamics of the edges. Interchannel impurity scattering does not influence 
strongly the tunneling exponents as they are found to be irrelevant processes at low energies. Exponents of 
backscattering at the point contact are unaffected by interchannel Coulomb interactions since all channels have 
the same chirality of propagation. 

PACS numbers: 73.43.Jn, 73.43.Cd, 73.43.Lp, 71.10.Pm 



I. INTRODUCTION 

An electron gas in two spatial dimensions (2DEG) in the 
presence of a strong magnetic field at particular values of 
the filling fraction (ratio of the number of electrons to flux 
quanta penetrating the sample) forms a charge gap in the 
bulk and gapless excitations at the edges creating there one- 
dimensional (ID) chiral metallic states. These states have at- 
tracted a great amount of interest as an arena for studying the 
exotic properties of the Luttinger liquid (LL) model'. 

Even though remarkable theoretical progress has been 
achieved in understanding the properties of Luttinger 
liquids, ^'^ in the experimental front open problems remain."*'^ 
For instance, questions have been raised from transport ex- 
periments, across a gate-created constriction in quantum Hall 
(QH) systems, by Roddaro et al.^. The measurements in 
Ref. la observe, in QH samples at integer bulk filling, a metal- 
insulator transition induced by varying the voltage of the 
metallic gates (V^) that define the constriction. At low |V^| 
the system is metallic (backscattering is suppressed) across 
the point contact and becomes insulating (backscattering is 
enhanced) as \Vg \ is increased. At high |V^| the observed cur- 
rent across the constriction has a bias-voltage evolution that 
would be reminiscent of a bulk-filling fraction v — 1/3. 

These types of tunneling experiments intend to probe pre- 
dictions of the LL theory such as the power-law dependence 
of the interedge tunneling correlation functions and the related 
power-law suppression of the densities of states. In QH sys- 
tems points of enhanced backscattering are created by bring- 
ing counter-propagating edges in close proximity by either 
placing metallic split gates over the top of the 2D electron 
gas, like the ones illustrated schematically in Fig.[Tl* or by al- 
lowing for breaks in barriers like the ones obtained by cleaved 
edge overgrowth^ (in which the interedge spacing is on the 
order of a magnetic length enabling additional electron cor- 



relations along the barrier). The theoretical groundwork for 
single-point tunneling phenomena in Luttinger liquids was 
laid by Kane and Fisher** by using a renormalization-group 
scheme. At low energies a change in nature of transport across 
a backscattering point depending on the value of the Luttinger 
liquid parameter was predicted: perfect transmission for at- 
tractive interactions across the sample and a perfect reflection 
at the point contact for repulsive interactions. These predic- 
tions were later confirmed by Fendley, Ludwig and Saleur— 
using the Bethe-ansatz technique which established the exact 
flow between the Dirichlet, T oo, the strong-tunneling sta- 
ble fixed point, and the Neumann, weak-tunneling unstable 
fixed point. Here T is the tunneling amplitude. Contrary to 
the experiments of Ref . 6 the LL theory predicts a rather unin- 
teresting voltage-independent top-bottom tunneling conduc- 
tance (G ^ y2i'-2-) ^j. fijjjjjg J/ = 1, typical for Fermi liquids. 
It is the intention of this paper to examine the disagreement 
between theory and experiment. 

The low jV^g I behavior of the differential tunneUng conduc- 
tance at a constriction in relation with experiments'" and— 
was examined in'^. It was concluded that repulsions between 
electrons on opposite sides of the gates suppress the quasi- 
particle (QP) backscattering at the constriction — thus leading 
to the metallic behavior observed at low | | . The idea there 
was that edges on opposite sides of the Hall bar must be re- 
garded as counter-propagating channels of a LL formed along 
the gates, rather than as the left and right portions of the top or 
bottom chiral channels of the overall Hall bar. Repulsive inter- 
actions between the chiral channels of such a LL are expected 
to enhance the left-right electron tunneling. 

The reason for the occurrence of the opposite phenom- 
ena: the enhancement of backscattering at the constriction 
observed at high values of |V^|, in systems similar to those 
illustrated in Fig. [T] was analyzed in Ref. [H. In this paper, 
we explain in significant detail the results that were discussed 
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FIG. 1 : Schematic illustration of a quantum Hall bar with a split- 
gate constriction similar to those used in experiments of Roddaro et 
al.^ . The electrostatic top gates bring counter-propagating edges in 
close proximity to enhance interedge backscattering. More details of 
actual devices are given in Ref.S 



in Ref. (Tl and elaborate on other experimental implications 
in connection with the proposed model. In analogy to the 
low jV^gl case, one might be tempted to attribute the exper- 
imental observations at high \Vg\ to the possibility of sign- 
reversed interactions across the barrier due to the presence of 
image charges in the metallic gates. To study the nature of 
interactions between edges on opposite sides of the gates, we 
allow for different potentials on the respective 2DEGs, and 
study the change of position of one edge upon a change of 
potential on the other 2DEG. We find in fact that electron in- 
teredge interactions remain weakly repulsive at any values of 
these voltages. The other venue that we pursue in this paper is 
the investigation of the effects of the edge-reconstruction phe- 
nomena that take place in wide and smooth edges inevitably 
created at high gate voltages, and caused by the presence of 
magnetic fields in QH samples. In the absence of a magnetic 
field we find that the width of the QH edge increases linearly 
on |V^|, with only logarithmic |Vg| corrections, favoring the 
need for an electrostatic treatment of such smooth edges. The 
charge distribution along the edges in the presence of mag- 
netic fields changes. With the introduction of charge energy 
gaps at points of fractional filling, the edge breaks into a se- 
quence of compressible (CS) and incompressible (IS) stripes, 
the latter exhibiting the fractional QH effect, QHE, (FQHE). 
Quite generally such stripes form in any situation where at 
zero magnetic fields a charge-density gradient is present, as 
noted in''* '^ '*. The ISs are strongly electron-correlated states 
and exhibit insulating properties as opposed to CSs which be- 
have as metal stripes (of constant potential in the area they 
occupy). The charge density remains with a gradient in the 
compressible ones and the low energy physics of the edge is 
described by their fluctuations'^. Large part of this paper is 
dedicated to the calculation of positions and dimensions of 
the CSs and ISs. The values of |V^| of Ref.|^for the geometry 
of these systems turn out to be such that for bulk v — \ filling 
the ratios of ISs to the CSs are of the order sa 1 and k, 0.5, 
for the outermost and the next to the outermost ones, respec- 
tively. The presence of these wide insulating stripes enables 
the backscattering of fractionally (e/3) charged quasiparticles 
whose tunneling exponent would correspond to an effective 
bulk = 1/3 filling. The interchannel interactions do not 
influence this exponent. The latter has its physical origin in 
the known fact that interchannel interactions cannot break the 
universality of the Hall conductance in same-chirality mul- 



tichannel edges like the ones at higher integer Landau level 
filling Hall systems. 

This paper is organized as follows. In Sect. II we examine 
the functional dependence of the edge width and edge posi- 
tion as a function of the voltage of the gates. Subsequently, 
in Sect. Ill we examine the nature of interedge interactions in 
the presence of metallic gates. In Sect. IV we find the geomet- 
ric dimensions of the compressible and incompressible stripes 
that form at the edges of the 2DEGs in the presence of a mag- 
netic field. In Sect. IV we find the geometric dimensions of 
the compressible and incompressible stripes that form at the 
edges of the 2DEGs in the presence of a magnetic field. In 
Sect. V we discuss implications of the reconstructed edge in 
transport and finite-frequency noise experiments. We close 
the paper with our conclusions. 



II. STRUCTURE OF EDGES IN SPLIT HALL BAR 
SYSTEMS IN THE ABSENCE OF MAGNETIC FIELDS 

In this section based purely on electrostatic arguments we 
study the charge distribution and edge structure and edge po- 
sition of a 2D electron gas in the presence of a confining po- 
tential provided by electrostatic split gates. We first examine 
the case of the absence of a magnetic field. The geometry of 
the model is represented schematically in Fig. [T] The 2DEG 
resides on the interface of two semiconductors, the xy plane 
(z = 0) of Fig. [U Two metallic gates are deposited over the 
top a distance w 10 /s from the 2DEGs' plane along the y axis 
of Fig. [T] Ib = \hc/Be\'^/'^ is the magnetic length, « 100 A, 
for a magnetic field B — 6.5 T that will be introduced later. 
The electron charge here is — e and c is the light constant. We 
assume here nevertheless that the gates and the 2DEGs are 
coplanar The gates, of lateral extension x E {—a, a), are 
used to create the constriction to enhance the interedge elec- 
tron tunneling amplitude. As the gate voltage is increased, the 
electrons under the gate are repelled leaving behind positively 
charged dielectric stripes parallel to the y axis of Fig. [T] ex- 
tended in the intervals a < |a::| < 6. In equilibrium, there 
is no net electric field on the 2DEGs' plane, and the 2DEGs 
can be treated as metallic components with constant potential. 
The equilibrium position will be reached at the points where 
the horizontal component of the electric field created by the 
charged dielectric stripes and the one created by the metallic 
components (the 2DEGs and the gate) sum to zero. At these 
points we require 

lim (Ef-{x,z = 0) + Ef-{x,z = 0)]^0 . (1) 

Each of those fields is singular at the sharp edges but in equi- 
librium they cancel exactly. 

Semiconductor structures in experiments'' were fabricated 
from AlGaAs/GaAs single heterojunctions of charge density 
no « lO'^'^ cm~^ and dielectric constant e = 12.6. In the 
depleted region the charge density is defined by the positive 
background charge of the donor semiconductor, cuq. In the 
2DEGs the electron charge density increases gradually away 
from the boundaries (located at ±6), reaching values of —eno 
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FIG. 2: Schematic illustration of split gates on the ?/ = plane at 
high gate voltages. For high values of | Vg| the metallic gate and the 
electron gases can be assumed to be coplanar. The potential of the 
2DEGs on the sides of the gate are taken to be equal to Vi. The 
positively charged stripes of the host semiconductor extend in the 
depleted regions a < \x\ < b. 



in the bulk. We are interested here in studying the case of 
high gate voltage, which justifies the assumption of putting 
the 2DEGs, the positively charged dielectric stripes, and the 
metallic gates all in the same plane, z — of Fig.|2] 

We assume here that the system is translational invariant 
along the gates. The problem then becomes effectively two- 
dimensional. A knowledge of the electric-field configuration 
in space leads to the knowledge of the electric charge distri- 
bution in space and thus to the knowledge of the structure of 
the edge which is the purpose of this section. We proceed 
in the following with the calculation separately of the electric 
fields created by the metallic components and the electric field 
created by the charged dielectric stripes. The potential on the 
plane z = fulfills the boundary conditions 



(X, z 



0) 



z) 



dz 




a < \x\ < b 



(2) 



(3) 



The latter is justified by the fact that the dielectric medium 
has a large dielectric constant of e = 12.6 [and 47r/(e + 1) « 
47r/e]. 

We solve in the following the Laplace equation in the half- 
plane z < using the boundary conditions dU and 
We separate the solution into a sum (j){x,z) = (j>'^^-{x,z) + 
(j)^*^'^' {x, z) of harmonic functions-!^ that satisfy separately the 
conditions 



0"'-(x,z = O) 
d(t>''^-{x,z) 



dz 



= 



\x\ < a 
\x\ > b 



a < Ixl < 6 



(4) 



(5) 



for the metallic components. The latter equation reflects the 
fact that the z component of the electric field in the depleted 
regions a < |x| < 6, for symmetry reasons, is zero. On the 
other hand, for the potential created by the dielectric stripes, 

(/>^'''- (x, z), we have 



dz 




t{x) 



\x\>b , 
a < Ixl < 6 



(6) 



(7) 



In these equations t{x) is the charge of the dielectric stripes, 

Aneno f : |x| < a or \x\ > b 



r(x) 



1 



a <\x\ < b 



and A{a, b) is a constant specific to our calculations and 
will be defined below. The boundary conditions fulfilled by 
(/i^'' {x, z) consist in the knowledge of the vertical component 
of the electric field in the depleted region and in the fact that 
the potential in the 2DEGs and gate regions is constant. 

To determine (lf*'''-{x, z) we generalize here an idea given 
by Glazman and Larkin"* by representing the potential 
(f)^*^'^- {x, z) as the imaginary part of an analytic function F{(), 
where ( = x + iz. By imposing the boundary conditions (|6]), 
(|7]i on the derivative dF/dC with an appropriate choice of an 
auxiliary function D{C,), Eq. (IA4b . we can relate it to another 
function, f{() = iD{(^){dF/dQ, whose imaginary part we 
would know everywhere along the x axis. By taking the op- 
posite steps (see Appendix [All we determine 



<j)'^'-{x,z) ^TbRe^E{k')- E{i;,k') 

-i±^[K(fc')- F(^,fc')]+i(^-fc 



(8) 



where ip — arcsin[ 



\J^J^y^ J"*"^^ and T — Aneno/e (more 
details are given in Appendix [All. Here and in the follow- 
ing K = K{k), E = E{k), K' = K{k'), E' = E{k') 
are the complete (complementary complete) Jacobi elliptic 
integrals of the 1*^' and 2"'^ kind of modulus k = a/b and 
k' — (1 — fc^)^/^,-'^ respectively. The potential constant 
A{a, b) introduced in (lU,® is given in ( |A9t . 

The potential (jf^- {x, z) created by the metallic components 
with the boundary conditions ^ and (|5]l can be determined 
by using 7r/2 folding conformal transformations at the points 
±a, ±6, mapping the upper half-space of Fig. |2]to a rectan- 
gle. To fulfill corresponding in the transformed space to a 
vanishing horizontal component of electric field on the boxes' 
vertical sides, the space inside the box has to be part of a larger 
space, obtained by successive iterative reflections of the box 
around its vertical sides, reaching eventually to an infinitely 
long capacitor whose plates are at potentials (lower plate) 
and Vi — A (top plate). The potential resulting from such a 
transformatior>2£ is obtained from 



Vi-A-V^ 



K{y/l-ayb^) 



c 



(9) 



with 0°' = Im$. The electric field created by the metal- 
lic components is obtained from E°^- — —lm[d^/d(]\z=o, 
whereas the one of the dielectric stripes is obtained from 
Ef - = -Im[di^/dC]U=o, cf. (|A8ll- Both of these fields ex- 
hibit singularities at the edge boundaries. To reach mechanical 
equilibrium at the boundary positions, one requires that 



Vi - K = 



Airenob 



K{y/l-a?lb^) - E{y/l-a?/b^) 



(10) 
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This equation^' defines the positions ±6 of the edges of the 
2DEGs. 

For small values of the dimensionless constant k :— {Vi — 
Vg)/Ta, with the expansion series for small values of the ar- 
gument in complete elliptic integrals,'^ one gets the linear re- 
lation b/a = 1 + 2k/tt. For the limit of large k, for which 
b — a ^ a, we make use of the logarithmic series'^ and get 
the asymptotic large |V^| dependence of b: 

ii^iVi~Vg)/Ta. (11) 



ln(4K/e) - lnln(4K/e) 

(Here the Euler number e is not to be confused with the ele- 
mentary charge e.) The distance of the edge position increases 
linearly on Vi — Vg with only logarithmic corrections. A plot 
of the dependence of the edge position as a function of the po- 
tential difference | Vi — Ig | obtained by the numerical solution 
of the transcendental equation ( fTOl i is shown in Appendix IB II 
(Fig. [IS}. 

The charge density distribution on the 2DEGs' regions, the 
z = plane at > 6, can be found similarly, by the knowl- 
edge of the z-component electric field. In the 2 = plane the 
stripes give a contribution of E^*''-'- = — Re [dF/dC] I^^q- 
derivative of F, with the appropriate behavior at a; — > ±oo is 
given by dF/d( = — irC''' (C) + t, where from Eq. ( IA8b 
we have G"*'' (C) = [C^ - (a^ + b^)/2]/D{C). In the regions 
> 6 we find 

N{x) = no{G'''-{x) - 1) - (^1 -A-Vg)G'''-{x), (12) 

which is the total charge density eN{x) — ~en{x) + eno, 
a superposition of the positively charged dielectric back- 
ground of density eriQ and the 2DEGs in > on top — 
carrying charge density —en{x). In Eq. ( fT2b the short notation 
G"=i (C) ^ b/[K'D{C)] is used. The profile of the edge of the 
electron liquid follows to be 



'n{x) = TT-o 



1/2 



\x\ > b 



(13) 



At large distances from the edge the charge density reaches the 
expected bulk value of no. The width of the edges increases 
as (6 — a), becoming smoother at higher \Vg\. Numerical es- 
timates for the edge position in devices used in Ref. 6, for 
- K = 1 V and a = 2500 A is 6 = 4964 A. The distance 
of the gate to the 2DEG is w 1000 A. This justifies our as- 
sumption of taking the 2DEG and the gate as being coplanar 



of nearby metallic gates is hard to estimate. The geometry ex- 
amined here, however, is such that it allows for exact calcula- 
tions. Intuitively one would expect that the presence of image 
charges in the gate, halfway the between edges, to weaken the 
interedge Coulomb repulsions or even reverse their sign. 

To study these effects we allow for different voltages on the 
2DEG subsystems. In the following we calculate the positions 
of the edges as a function of the voltages and look at the ra- 
tios of the change of the position of one edge upon the change 
of voltage on each of the 2DEGs. Due to the voltage differ- 
ence between the 2DEGs, their distance from the metallic gate 
should be different. We assume in general the edge boundaries 
be positioned at b and b whereas the gate be extended within 
the interval a < x < a, with a — b ^ b — a. 



A. Field of Asymmetric Dielectric Stripes 

Here we calculate the electric field and the potential config- 
uration created by the positively charged asymmetric stripes 
and its constant potential surroundings. The dielectric stripes 
extend in the intervals (b, a) and (a, b) as shown in Fig.[3](up- 
per left). The boundary conditions for the dielectric stripes 
now have to be modified to 




a < X < a 
X <b , 
X > b , 



(14) 



t{x) , b < X < a , a < X < b , (15) 



where t{x) — Aneno/e and A, B are constant potentials that 
will be defined below. 

We impose the boundary conditions (O, ( fTSI l on dF/d(l, 
where 0*'*'' (a;,2) — ImF(C) and C, — x -\- iz. We intro- 
duce as before the auxiliary function D{C,) (Ref. l23l) to build 
/(C) = \D{C,)dF / dC, such that /(C) takes imaginary values 
everywhere on the x axis. Using the Schwarz equality ( IA5I ), 
we get in analogv to_(lA6Ti: /(C) ^ t[iD{C) + T{Cb,a,a,b)]. 
The function J-{C,, b, a, a, b) again must be chosen with the 
correct behavior at infinity, /(C) as |C| — > oo, resulting 
in !F{C, b, a, a, b) given in^^ for constant t. The knowledge 
of !F{C, b, a, a, b) immediately leads to the knowledge of the 
electric field in the depleted regions, z = 0, b < x < d, 
a < X < b. 



III. NATURE OF INTERACTIONS 

In this section we study the influence of the metallic gates 
on the electron interactions between edges on opposite sides 
of the gates. Previous work'^ has found that QP or electron 
backscattering in the constriction can be enhanced or sup- 
pressed depending on the sign of the interedge interactions in 
the constriction. Attractive interactions, for instance, would 
effectively lower the filling fraction in the Hall bar and en- 
hance the backscattering of the emerging QP at the point con- 
tact. In general cases the influence on the electron interactions 



J^{x, b, a, a, b) 



(16) 



The electric field created by the stripes, as one would expect, 
is singular at every edge (including the edge boundaries b and 
feofthe 2DEGs). 

The potential created by the dielectric stripes, c/f^'^- (C), after 
integration of /(C)/-D(C) (Ref. l24l) is found to be given by 
the expression ( lAl lb . The constants A and B are defined as 
potential differences A = 0''*'' (6, 0) - (a, 0) and B = 
i/i^'' ■ {b, 0) — (a, 0). From ( lAllb the expression for B turns 
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FIG. 3: Upper left: original potential configuration in space — metallic gate and 2DEG potential. Since Vi 7^ V2 the distance of the 2DEGs 
from the gate is different. Lower left: after inversion around a point si the 2DEGs' distances from the gate are equal. Right: forming a box 
configuration (shaded) by mapping the upper C,i half-plane, which in turn maps the lower half-plane of the original plane Application 
of the Schwarz mirroring principle about the vertical borders of the shaded box leads to a periodic structure along the X2 axis. 



out to be 



B = -SI 



b~b 



a — a 
b-b 



K{k')--S+E{k'), (17) 



where the modulus of the elHptic functions is [trans- 
formed from k of (lAlll i to] k = Sj^/S'^ with = 

[{b — a){a — 6)] ^^'^ ± [{b — a) (a — 6)] The expression 
for B is invariant under b ^ b, a a, implying that A ^ B. 
This constant represents the value of the potential at large dis- 
tances from the gates. 



B. Field of Gates and 2D Electron Gases 



(lower left). Finally we choose s such that the distances be- 
tween the 2DEG edges and the metallic gate in space be 
equal, 1(a) — I{b) — 1(b) — 1(a). This condition is resolved 
at two points, si_2r^ from which only si is used below. The 
transformation in (l20l l inverts the space bringing the points 
±ioo to =FiO and ±00 to ±0. In Ci space the gate lies in the 
outer regions and the dielectric stripes around the origin en- 
closing the 2DEGs. The 2DEG subsystems of potentials Vi 
and V2 are located, respectively, between /(oo) and 7(5) for 
Vi, and between 1(b) and /(oo) for ¥2- 

Second, we perform another conformal transformation 
mapping the new configuration to a box [tt/2 folding trans- 
formations at the points I(ai), i = 1, . . . , 4]. This brings us 
in the C2 plane. 



As in the previous sections, we assume that the metallic 
gates and the 2DEGs are coplanar. The boundary conditions 
for the potential of the metallic components now have to be 
modified to 



(t)°'-(x,z^ 0) = <( Vi=Vi-A 

V2=V2-B 

i'='-(a;,z) 



a < X < a 
X <b , 
X > b , 



(18) 



dz 



: b<x<a,a<x<b. (19) 



The potential difference of the 2DEGs renders the tt/2 fold- 
ing transformations, used previously, not directly applicable. 
We circumvent this limitation and determine the field config- 
uration in space by employing the following steps. 

First, make an inversion transformation around a point s on 
the X axis under which the points ai, i = 1, . . . ,4 correspond- 
ing to 6, . . . , 6 of the ( = X + iz plane map to points /(a^) as 
follows. 



C2 



Jib)' 1(a)) 



\Iia)J 



(21) 



where we have selected the multiplicative and the additive 
constants be 1 and K, respectively, as shown also in the right 
panel of Fig. [3] The modulus of the Jacobi elliptic integral 
K is 1(b) 1 1(a), which equals k — S/^ / of the Jacobi inte- 
grals contained in the expression for B of the previous section. 

To fulfill the boundary conditions ( fT9l ) of vanishing electric 
field along z in the depleted regions in C, space, the shaded 
box of Fig. [3] (right) in (2 space has to be part of the periodic 
configuration shown on the same figure, in which the electric 
field vanishes on the vertical borders of the shaded box. 

The expression for the potential field in the new space can 
be calculated analytically to the close form (see Appendix 



2K 



K' 



Ci 



i?2 



p 



Ha,) 



i?2 



P 



(20) 



R, the radius of inversion, introduces an overall length scale. 
The potential field is scale invariant and independent of R, 
however P merely shifts the origin of the d plane, which we 
take to be the middle point between 1(b) and 1(b), see Fig. [3] 



+ — ^Re<iln 



h{-ik(C2 + c)., 



e 2K 



^1 



{jkiC2-C),e- 



(22) 



Here -di is an elliptic theta function'^ and C is the reflection 
of the ±00 points of ( space in (^2 space, C = C2(/(±oo))- 
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This field configuration can now be mapped onto the original 
plane. 

Now we have the necessary ingredients to proceed on de- 
termining the functional dependence of the 2DEGs' edge po- 
sitions b, b on Vi and V2 ■ 



C. Edge-state Position as a Function of Potentials and the 
Influence of Metallic Gate on Interedge Interactions 

Since the edges of the 2DEGs are mobile, their equilibrium 
position is reached at points where the x component of the 
electric field in the ^ = plane of the overall system vanishes. 
Once again we use the mechanical equilibrium equation sim- 
ilar to Eq. ([B: 

E,{x,z^O)l^-^^0 , E,ix,z^O)l^^^O ,(23) 

and try to turn these conditions into a system of equations 
for the edge boundaries b, b. The x-component electric field 
at a; ^ 6 and x ^ b in the z = plane of both, the di- 
electric stripes, as well as the metallic components in princi- 
ple are known. The calculations in this section consist pri- 
marily of a simplification of the expressions for the electric 
field of the metallic components at z — 0. We seek a refor- 
mulation of from the elliptic di functions, of modulus 
exp{—TTK' /2K}, to elliptic integrals of the same modulus 
as the ones of the expression for the potential B of Eq. ( fTTI i 
that is present in the potentials of the capacitor plates. For the 



sake of clarity we give some details here, but leave the rest to 
Appendix lB 21 

First, for the electric field of the metallic components the re- 
lationship between E^^- and Ey]^ at the points ^2 = 2nK + iy2 
[corresponding to points in the intervals x S (b, a) and x E 
(a, b) and z = in C space] is needed. Using the equalities 
( I20I ), tni one can show that at the points (2 = 2nK + \y2, 
where the field does not have a horizontal component in (,2 
space, the Jacobian of the transformation C ^ C2 takes a sim- 
ple form, and hence 



5C2 dCi 



■ (24) 



X2=2nK 



For the coordinate transformation (|24] |. with the explicit form 
of si from Ref.Hil one obtains: 



.9C2 5Ci 



J(2^0) = -i^^ ^ 



1 Si 



2DiC) 



(25) 



The Jacobian of the transformation to the original space pro- 
vides the expected square-root singularities for the x compo- 
nent of the electric field at the edge boundary points b and b. 
The sign of the electric field is determined by the factor -D(C) 
along the line ( = x — iO. 

Second, the vertical component of the electric field in (^2 
space is given by 



J 



£^V2 (C2, C2) 



gr-(C2,C2) 

dy2 



-Vo 



1 V1-V2 
K' 4K 



Re 



,^1 (4^(C2 + C),e-^ 



(26) 



where the prime in denotes the derivative with respect to 
the first argument of andVb = Vg_-V2-{Vi-V2)C/{2K) 
was introduced. 

The conditions ^ involve the fields E^^ at C2 = 0, 2K, 

dMJ, (Hg. With dmOb to dsni l and by using Vi instead of V, 
and Vq = Vq — B one gets for their symmetric and antisym- 
metric combinations. 



b + b 



El\{2K) - E';}^{Q) 



y2 ' 



(Vi ~ V2) (b~b) 
_2li , (27) 

Etim+EtM = - , (28) 



K' 



where the length I is given by Ref. [23 On the other hand, 
the field created by the charged dielectric stripes is obtained 
from (fTSl i. Inserting these field expressions and (fTTI i into the 
difference and sum of fields obtained by multiplying (|23l l by 

D{x), we get 



Vi~V2 



TTT 



E' , 5+^0 , V1-V2 



K' 



(29) 



I (30) 



with (5"+)" < (5 - 6)2 and E' < K' . Equations ^ and 
( l30l ) define a complete system from which the dependence 
of the edge positions b and b of the 2DEGs on the voltages 
of their respective potentials can be found. We have solved 
these equations self-consistently numerically and have plot- 
ted results for 6(Vi, V2) for fixed Vg in Fig. |4] The plot 
for 6(Vi,V2) is equivalent to V2) under the exchange 
Vi ^ V2, X <-> —X. In the limit when the voltages on the 
2DEG subsystems are equal, Vi = V2, we have a — —a, 
b = —b and get S'^ = 2b. The equations for the edge po- 
sitions in this limit simplify, yielding an equation consistent 
with the earlier result [Eq. ( flOb l that was derived slightly dif- 
ferently in Sect. HIl 

The quantity of importance in our calculations here is the 
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G 



FIG. 4: Position of tiie edge, 6, on tiie right of tiie metallic gate 
as a function of Vi — Vg and V2 — Vg. At negative gate voltage 
h increases past the edge of the metallic gate. The dependence for 
fixed Vi is linear on V2 — Vg. For Vi equal to the distance of 
the edge positions from the gate is linear in agreement with the one 
shown in Fig.[T6l 




FIG. 5: dvib/ dv^b, the ratio of the differential dependences of the 
edge position on the potentials of the 2DEGs on the sides of the gate. 
By increasing Vi of the 2DEG subsystem on the left of the gate of 
Fig.|2l the edge on the right-hand side approaches the gate, just as the 
other one on the left goes away from it. This signifies that interedge 
interactions in the presence of the metallic gate remain repulsive in 
nature. 



sign of dvib/dv^b. By increasing Vi, as can be seen in Fig.|4] 
\h\ increases, as does h upon the increase of V2. The graph 
in Fig. |5] shows that upon the increase of Vi the 2DEG edge 
6, on the other side of the gate, approaches the gate as the 
other edge h goes away from the gate. This signifies that 
interedge interactions remain repulsive in nature at any val- 
ues of Vi and V2. At low values of Vi — V^, i — 1, 2, in- 
stead dvib/dv^b ^ 0, signifying the tendency for the com- 
plete screening of interedge interactions from the gate. We 
conclude that the explanation of enhancement of top-bottom 
backscattering (the insulating behavior across the constric- 
tion) observed in the experiments Ref. 6 has to be related with 
other phenomena taking place at high gate voltages. We elab- 
orate on this in the following section. 



IV. STRUCTURE OF THE EDGE IN THE PRESENCE OF 
MAGNETIC FIELDS 



n(x) 



n(x) 



n(x) 



-b b 



FIG. 6: Structure of the edge of a 2DEG surrounding the split gate, 
in the presence of a strong perpendicular magnetic field in QH sys- 
tems with constrictions. Shaded areas notify current-carrying CSs, 
separated by insulating ISs. For half a gate width of a = 5000 A 
the geometrical dimensions of the stripes are 180 A for the first CS, 
180 A for the 1/3 IS, 1040 A for the second CS, and 377 A for the 
2/3 IS, respectively. 



ofRefs. [l4lfl5lfl6i andlH In the edge of a bulk integer = 1, 
around points of / — p/ (2mp ±1) filling, for p and m in- 
tegers, a hierarchy of FQH energy gaps exists. The filling 
fraction as usual is defined hy ly = no/27rZ^ as the num- 
ber of electrons per flux quantum. With the introduction of 
the gaps, the potential in the 2DEGs is not constant through- 
out as in Sect. |ll] As we move in the x direction away from 
the edge boundaries, the potential varies in a discontinuous 
manner, jumping in value by Vf each time we cross the lo- 
cation of a FQHE gap corresponding to filling /. The sys- 
tem therefore will tend to lower its energy if some electrons 
are relocated from higher to lower energies forming this way 
fractional-filling ISs around points of f ~ p/{2nip± 1), each 
of which exhibits the FQHE. The energy profile across the 
edge varies in a step-like manner, being constant in the com- 
pressible metallic regions and increasing across the incom- 
pressible insulating stripes. Therefore application of a strong 
magnetic field will in principle break the edge of the 2DEG 
confined by a smooth external potential into a sequence of 
CSs and ISs. The widest ISs will be the ones with the largest 
gaps [see also Eq. ( l45b below] and these will be the ones real- 
ized in practice. In the case of integer fillings the role of the 
FQH energy gaps will be played by the much larger Landau 
level splitting. As a result the integer edges will lead to wider 
and more stable incompressible stripes. The positions of the 
incompressible stripes correspond to fractional fillings and in 
general follow equipotential lines parallel to the external gate 
(Fig.©. 

The largest charge gaps in the FQHE are known to occur 
at / — 1/3 and 2/3 fractions^^ and in the following we cal- 
culate their widths. In general the position of a fractionally 
charged stripe of filling fraction / can be found by substitut- 
ing n{xf±) = /no in Eq. ( fT3l ) and solving it for Xf±. Such 
substitution leads to the Xf± positions 



The effect of strong perpendicular magnetic fields on the 
structure of the edge of a 2DEG is known from earlier work 



Xf± 



= ±6 



1 - Pf 



1/2 



(31) 
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Based on the estimate of Sect.|II]for b = 1.986 a at Vi 
0.674 V and B 6.5 Tesla, a = 25 we find 



Xl/3± 



±2.08 a 



X2/3± 



±2.51 a 



(32) 



In this section from the knowledge of V/'s we calculate the 
widths of the sequence of CSs and ISs on the edge. It is as- 
sumed that the metallic gates and the 2DEGs are coplanar. 
The electrostatics problem in the presence of a magnetic field 



.15 



(j)'{x,z^ 0) 



d(j)'{x, z) 



: |a;| < a 
Vi : & < |a;| < a' 
Vi + Vf : \x\ > b' 



dz 



: a<\x\<b, 
'-{hq ~ Uf) : a' < \x\ < b' 



(33) 



(34) 



In ( [34] i a factor of two is missing when compared with (O. 
This is due to the fact that the ISs are narrow when com- 
pared with the 2DEG distance from the surface of the semi- 
conductor, and the electric field surrounding the ISs can be 
thought of as concentrated within the semiconductor. In (|33] | 
Vf is the voltage difference on the sides of the IS of filling 
fraction / that arises due to a discontinuity on the value of 
the chemical potential of elctrons across the IS. Vf is given 
by Vf = A^f/e, where for integer filling IS, A/i — fiUc, 
whereas for fractional fillings / = p/q, where q is an odd 
integer, it is given by Afif — o cfe^/eZe. The constants Cf 
are Ci/3 = 02/3 = 0.03 (Ref. 1271 and references therein). The 
quantity Cfe'^/els, the charge gap, is the energy necessary to 
create a pair of (particle-hole) quasiparticles of charge e/q. 
The biggest charge gaps are known to occur at / = 1/3 and 
2/3 fillings. In the following, therefore, due to the smallness 
of the charge gaps, the ISs at other fractions will not be con- 
sidered. For / = 1/3 and / = 2/3 we have = 1.026 mV. 
As in Ref. iTj the solution of the Laplace equation in the half- 
plane z < can be found as a sum 0' = 0'°'- + + (pQ 
of harmonic functions 0'^'"' and the zero-magnetic-field 
solution 00- We model the boundary conditions for (/)'°'' as 
the ones of the metallic 2DEGs of Sect. and those of (Z)"**' - 
as the ones belonging to the corresponding charged dielectric 
stripes. The boundary conditions for then read 



(/)"='-(a;,z = 0) = 



d(l)"''ix,z) 



dz 







: |a;|<a', 
VI : |x|>6', 



a' < \x\ < b' 



(35) 



(36) 



where V{ — Vf — A{a' , b'). In analogy to the boundary con- 
ditions ©, © for (lf^^-{x, z) the potential (j)'^'^'' ■ {x , z = 0) of 
the charged stripes formed at the edge profile fulfills 



■{x,z^{)) = 




\x\ < a' , 
\x\>b', 



(37) 



■{x,z) 



dz 



27re 
= [n{x) 



f], a' < \x\ < b' . (38) 




0.5 1 1.5 2 
V,= 1.026 mV > Vi-VgCV) 



Aa,/i (A) 




V,/i 



(Vi-Vg)/x (A) 



FIG. 7: Plots of width of incompressible stripes of 1/3 filling Aai/3 
in A as a function of Vi — Vg and Vf for the second-order approxi- 
mation and the values Vf (i/3) V/, i = 1, . . . , 5, from bottom to 
top are shown (in Ref. 6, Vi ~ Vg ^ 0.674 V and Vf = 1.026 mV). 
The gate width is taken to be 2a with a = 251b- 



The only difference with the previous case is the x dependence 
of the charge density of the stripes extended in a' < |x| < 6' 
on the sides of Xf. Qualitatively, a presence of the potential 
difference Vf on the sides of Xf will tend to push away the 
sides, and the remaining stripes will provide the traction for 
the electrostatic equilibrium at the points a' and 5'. The fact 
that the stripes are now charged both positively on the left of 
X f and negatively on the right of Xf will result in a reduction 
of the traction at small Vf as compared to the case of uniform 
positively charged stripes provided by the dielectric of Sect. 
nil Al leading to a rapid opening of ISs at small Vf. Neverthe- 
less, as Vf increases the positive part of the charged stripe will 
become strong, leading to an eventual slowed-down increase 
of the width of the ISs as a function of Vf at higher values of 
Vf- 

Since the values of the charge gaps are relatively small 
when compared with |Vi — V^|, V// | Vi — « 10"'^ in 
Ref. m the charge density in the region of ISs, N'{x) — 
n{x) — rif, can be expanded around the points Xf±. The 
even function n(x) is approximated by one even polynomial 
J2iLo nf^^x^'^ matching n{x) up to the iVth derivative at x f±. 
In the following the notation t{x) = {2'Ke / e)N {x) is used. In 
electrostatic equilibrium the positions of the IS edges are ob- 
tained by requiring the electric field created by the metallic 
components extended outside the intervals a' < < 6' and 
the charged stripes, be zero at all points ±a', ±6'. 

In principle the calculations for the electric fields E'^'^-{x) 
and E^^'^- {x) satisfying the boundary conditions ([35ll-(l38]l for 



i/ol. 



and are similar to the ones of Sect. [Ill The 
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only difference arises from the fact that the charge density of 
the stripes has now an a;-dependent profile. In the follow- 
ing we describe briefly the calculations for i?*^*' (x), leaving 
details of the calculation in Appendix IA2I By representing 
the potential (p^*"''- as the imaginary part of an analytic func- 
tion F, and using the boundary conditions for its derivative, 
we relate this to the function /, /(C) ^ iD{Q{dF/dQ, 
( = X + iz, which in turn after using Eq. iA5\ takes the form 
/{() — T[iD{() + ^(C)], where the differences with the pre- 
vious case come primarily in the form of !F{C,). will de- 
pend on the form of the D{() as well as on dF/d^, i.e., on the 
function t{x) = J^n '^nX^'' of (l38T l. where n will be the order 
of our approximation. The general form of T{C,) — X]„ 
results form the knowledge of 

[n/2+l] ^ 
fc=0 

(39) 

where Cj, ^ (t) is a Gegenbauer polynomial, and [•] represents 
the integer part of a real number. is found straightfor- 
wardly by integrating f{Q/D{(). 

In the following we will calculate the cases of n = 2 and 
n = 4 of the approximation. 



1. Second-Order Expansion 

Here we take the lowest order of approximation for n{x) of 
(O, i.e., N = 1, 



(40) 



sity of the stripes extended in the intervals a! < \x\ < b' 
will change the results slightly as compared with the ones of 
Sect, ini The result in this case for the previously introduced 
F function, whose imaginary part gives the potential, is ob- 
tained by the integration of ( IA8I ) and given in ( lAl II ). For A in 
Eq. (l37b we get 



A(i)(a', b') = 4^>b' (i(l + k^)K' - E') 



(41) 



with r^^^ = -x/t72, t^^^ = t' /2xf, t' = 27ren|^/e, and 
k — a! jb' . As in Sect. |II]the value of the potential constant 
yl(a',fe') is specific for our calculations. Requiring for the 
total value of the electric field at both edges of the stripe to 
vanish, we obtain two equations, one of which reads 



111,12 ri\ 

T (b — a ) 



an 

4 



36'- 



16xf 



Vf-A(^\a',b2^, 



K' 



(42) 



and for the other one a' and b' are interchanged on the left- 
hand side. Their numerical solution would lead to the stripe 
width Aa j = b' — a'. From the difference of ( |42] ) with its left- 
hand side (a' ^ b') interchanged a constriction for the edge 



boundaries of the ISs is gained. 



2xj 



(43) 



Xf is away from the middle point between a' and b' shifted 
toward b', a feature that becomes more pronounced for wider 
ISs. Taking the sum of (l42l l with its left-hand side (a' ^ b') 
interchanged using also ^(i'(a',6') from gD we obtain as 
second equation for the position of the stripe boundaries a' 
and b': 



F2 



E' - k^K' 



3 1 



k'^y Vf 



(44) 



At the small values of Vf considered here the above equality 
therefore will be resolved at values A; ^ 1 (6' — a' <C a'). 
In this limit using expansions of E' and K' for small k' 

for fc we get fc w 1 — {2/ x f){2Vf /ttt'Y^'^ , which gives 

1/2 

Aai/3 « [4eV/-/7r^e9a;n(a;i/3)] . The stripes will be wider 
for bigger gaps and around points with smallest charge gradi- 
ent. Therefore the ones located further from the edge bound- 
ary (inner) ISs will be wider than the ones located closer to 
the boundary (outer ones). 

Because asymptotically Xf increases linearly with 6, n'f 
decreases accordingly, and so the term l/(T'a;j) from ( |44| ) 
decreases inverse linearly. Therefore its solution has to be 
taken in the limit k ^ 1 ib' — a' <^ a'). The widths of ISs, 
Attf = b' - a', turn out to be^° 



with n'f = dxn{x)\^^. The x dependence of the charge den- 



1 fVf 



Treno \tt (1 _ p)i In [A{Vi - Vg)/Tae\ 



,(45) 



proportional to (Vi — VgY^^, up to logarithmic corrections. 
The dependence of the IS widths is also proportional to 

1/2 

Vj: , in sharp contrast to the previous case, of the uniformly 
charged stripes. Such a behavior has its origin in the fact that 
the stripes have a charge profile in x direction. Plots of the 
dependence of Aa/ onVi—Vg and Vf, obtained by the nu- 
merical solution of ( l43T l and ( l44b . are shown in Figs. |7] and [8] 
(top). 

For the narrower outer ISs we find a width of w 165 A- 
180 A. The wider inner 2/3 channels, of size w 340 A- 
w 380A, are « 800-1000 A away. 

For the additional charge, induced by the magnetic field in 
forming the ISs at | a; | > 6, we find 



An = 



2x 



\D'{x)\ 
-\D'ix)\ 



: \x\ > b' 
: \x\ < a' 



where D'{x) is obtained by D{x) of Sect. |ll]by the substitu- 
tion a ^ a' and b ^ b' . A plot, on a magnified scale, of An 
along the system of two edges is shown in Fig. [8] (bottom). 
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FIG. 8: Top: three-dimensional plot for stripe widths Aa/ for 
/ = 1/3, measured in units of half of the gate width, a, taken usually 
in Ref 6 to be a = 2500 A-5000 A. Bottom: Schematic plot of the 
additional electron charge caused by the presence of a strong mag- 
netic field on the edge profile of a QH state wide enough to support 
ISs. a — 3a, b' = 5a are chosen as the positions of the boundaries 
of the 1/3 IS on the right of the gate in this illustration. The point 
a' < Xf < b', where n{x /) = /no, is defined in the text. 



2. Fourth-Order Expansion 

Here we require in the approximation also the second 
derivative of n{x) at a;/ to coincide, 



n{x) 



1 



1 



n'i ] x^ 



1 



(46) 



where 



d'^n{x)\xf The equations for the stripe 
edge positions will be given by a set of equations similar 
to ( |42] |. where the denote the coefficients of t{x) — 

J^iLnJj^^ ^'^^ (with N — 2 here) in Eq. dSST l according to 
Eq. (|46|: 



-(2) 



J2)a 



'2 + 3b''- 



+ 2a'H"' + 5b'^ 



Vf^A^'\a',b') , 
K' 



(47) 



where again the second equation follows by the exchange 
a' ^ b' on the left-hand side of Eq. (|47] ). Difference and sum 
of these equations are given in Appendix IA 21 Their numerical 
solution is shown in Figs.|9](top) and[TT| A solution exists for 
not too large values of Vf because of > and n'j < 0, 

(2) (2) 

hence > and ' < 0. This is because the fourth order 
of approximation of n{x) has an upper bound (as the original 
profile Eq. ( fTsT l, but presumably at slightly smaller values) in 




FIG. 9: Top: three-dimensional plot of the dependence of the 
incompressible-stripe width Aa/ for / — 1/3, formed at the edge 
of the 2DEGs in the presence of magnetic fields, on gate voltage 
Vi — Vg and charge gap Vf for the fourth order of approximation 
of n{x) around points Xf±. For large values of Vf no solution ex- 
ists. Bottom: Schematic illustration of fourth order for An. (a' = 3, 



b' = 5, coefficients of r(x): r^'^'' = 1.25, r^^' = -0.025.) 



.(2) 



n 








1 1 


b' 




:'1.6 


\Y8 1 i^^^i-— 



FIG. 10: Schematic illustration of second (gray) and fourth order 
(black) approximation for An around the right IS. a and b' are 
numerically determined for magnified values of n{x), Vf/r'. The 
fourth order of approximation leads to slightly wider ISs than the 
second-order one. 



contrast to the second order ( |40] |. which is unbounded at large 
values. Therefore, a solution for Eq. (l42b can be found for any 
positive Vf, whereas for Eq. ( |47l ) it cannot. 

In fourth order the additional charge An is given by 



An : 



^(2) 

^2 



1 

TO 



-fe'2 



\D'{x)\ : \x\ > b' 
-\D'{x)\ : \x\<a' 



with To = 27re/e, a plot of which on a magnified scale along 
the edge is shown in Fig. |9] (bottom). 

In the following we give comparisons of sizes of incom- 
pressible and compressible-stripe widths between calculations 
of the second and fourth order. These two forms of approxi- 
mation for Ati are plotted in Fig.fTOl 

The fourth-order calculations depart from the second-order 
ones only signifcantly when V//(Vi — V^) « 0.01, Fig. [TTI 
The experimental values of the ratio of these voltages are 
Vf /{Vi — Vg) ~ 0.0015, and so the fourth order will not give 
appreciable corrections as shown in the table below. 
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(l/-y)/(Xfl)-1.86 Vl{%a) 




0.5 1 1.5 2 
VAxa) = 0.002832 (V - V )l{xa) 

J '- o 



FIG. 11: Dependence of the width Aa/ of the incompressible stripe, 
of filling fraction / = 1/3, on the magnitude of the charge gap Vf 
(left), and on Vi — Vg (right) for second (triangles) and fourth order 
(squares) of approximation of the charge density around points xf±. 
Notice the square-root dependence of Aay for small Vf. At large 
gate voltages, Aa / also shows a square-root dependence on Vi — Vg. 
A solution in fourth order for large values of Vf compared to Vi — Vg 
does not exist, see also Fig.|9] 





a = 2500 A 
Vg = -0.674 V 


a = 5000 A 
Vg = -0.674 V 


v = 1 


6 = 4964 A 


6 = 7649 A 


/ = l/3 

2nd 

4**^ order 


Xi/^± = ±2.076 a 
Aai/3 = 164.2 A 
Aai/3 = 165.4 A 


xi/^j^ = ±1.544a 
Aai/3 = 179.4 A 
Aai/3 = 180.3 A 


/ = 2/3 

2nd Qj.jg^ 

4*^^ order 


X2/3.± = ±2.509 a 
Aa2/3 = 338.0 A 
Aa2/3 = 338.9 A 


X2/3± = ±1.847a 
Aa2/3 = 375.9 A 
Aa2/3 = 376.8 A 



The fourth-order approximation, however, is expected to give 
noticable contributions in the case of edges of higher Landau 
level fillings when V/ is one order bigger than the case of the 
integer filling discussed here. 

If the voltages V/ were on the order of Landau level split- 
ting (about 10 times larger than the current V/), which is 
the case of edges of higher integer filling fractions, then the 
IS widths in 2"*^ order of approximation for / = 1/3 and 
/ = 2/3hadasize of519.6Aand 1071 A, respectively. There 
the 4**^ order would lead to sizable effects, leading to 582.6 A 
and 1106 A, respectively, cf. also Fig.fTTl 

Because b grows almost linearly with a, the changes of the 
widths of the ISs are rather moderate as a is varied. 

In experiments'' the split gates are slightly out of the 2DEG 
plane (w 10 1 b off plane in z direction). As a consequence the 



gates provide a smoother confining potential than the copla- 
nar case examined here. This should be expected to lead to 
smoother edges and wider ISs. 

V. EXPERIMENTAL IMPLICATIONS 

In the following we discuss implications arising from, re- 
spectively, the composite multichannel edge structure sur- 
rounding the split-gate constriction at high | | in QH systems 
at integer bulk fillings, and the inhomogeneously interacting 
non-chiral LL system in which they turn to at low |V^|, in 
transport experiments. Their observation would offer a fur- 
ther proof of the consistency of the model that we propose to 
explain the nature of the metal-insulator transition observed 
in Ref.'el 

At low I I , measurements performed over a broad range of 
filling fractions 1/3 < < 1 (Refs. 6 10 11) consistently ob- 
serve a QP backscattering suppression across a constriction. 
At these voltages the separation of counter-propagating edges 
under the gate is small and consequences of interedge inter- 
actions correspondingly are expected to be most significant. 
The mechanism that we proposed therefore in Ref . [Tj and ap- 
plicable for experiments^ where v ^ \, attributes backscat- 
tering suppression to interactions near the point contact. This 
can be most easily understood by regarding the edges on op- 
posite sides of the line junction as counterpropagating chan- 
nels of a non-chiral LL defined on the constriction. These ID 
systems form an inhomogeneous LL with parameter K < 1 
in the junction region (enforcing relevancy of left-right elec- 
tron tunneling or suppression of top-bottom QP tunneling) 
and K ^ 1 in the remaining portions of the edges. These 
systems have properties that can be observed in experimental 
studies as we elaborate in Sect. IV Bl below. 

At high |V^| it is found in Sect. |III] that interactions across 
the line junction cannot lead to the opposite phenomena of en- 
hancement of QP backscattering. We explain below that in the 
reconstructed edge at high | Vg | relevant fractionally charged 
QP tunneling processes arise and an insulating behavior is to 
be expected. Our proposal for the explanation of experimental 
observations at high |V^| is captured by the schematic illustra- 
tion represented in Fig.[T2l 

We comment in the following also on the features one 
should expect from noise experiments on these systems at 
high and low |Vg|. 

A. High Gate Voltage 

Detailed calculations in Sect. |IV] found that for the geome- 
try of the experimental system and the values of |V^| that are 
used, the 1 /3 and 2/3 ISs are of widths 2 1 b and 4 Ib, respec- 
tively, separating CSs of sizes 21b and 10 Zs, respectively. 
These ISs are wide when compared with the width of electron 
wave functions {I b for the noninteracting case but less than I b 
for the interacting one) in lowest Landau level states. There- 
fore wave functions for electrons in guiding centers localized 
in front of each other on opposite sides of the ISs, have an 
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v=l 



FIG. 12: Schematic illustration of high | Vg \ splitting of a = 1 edge 
into channels of fractional filling in QH systems with constriction, 
similar to those employed in experiments of Ref.@. 



overlap of the order of e 1/3/ « ^ where Aa^ ^3 is the width 
of the narrower 1/3 IS, leading to small inter-compressible- 
channel tunneling amplitudes. Moreover, wide ISs enforce 
an effective statistics for the fluctuating CS channels similar 
to the one characterizing the v — 1/3 edges. In this pa- 
per, therefore, we propose the model depicted schematically 
in Fig. [12] as the one that develops in the experiments at high 
|V^| and that explains the insulating behavior across the con- 
striction. In the absence of interchannel Coulomb interactions, 
the LL formed between ly — 1/3 channels at the line junction 
anticrossing the Hall bar is equivalent to a LL along the junc- 
tion at filling l/i^,'^'"" a reflection of quasiparticle-electron 
duality in tunneling processes. The flow of G from high to 
low V can be described as follows: Starting with an open con- 
striction (Fig.fTSTi. slowly turning off V, the system flows from 
relevant quasiparticle backscattering, suppressing the source- 
to-drain conductance according to G — ve^ /h ~ —V"^"^"^, to 
a closed constriction with irrelevant electron tunneling with 
open the constriction at small V . In 
the absence of interactions the system is insulating at low bias 
with I-V characteristics of a QH system of bulk v — 1/3. 
The I-V evolution being reminiscent of a bulk filling 1/3 was 
indeed one of the puzzles of the experiments^^ 

The chiral channels constituting the integer edge are in 
close proximity, nevertheless, and interactions cannot be ne- 
glected. We show now that Coulomb interactions between 
channels do not influence the tunneling exponents. In the fol- 
lowing we adopt a model of coupled chiral channels, each de- 
scribed by a single chiral LL, an assumption best fulfilled by 
the outermost sharper channels. The excitation spectrum of 
this model is composed of two slow neutral modes and a sin- 
gle faster charged mode, each extended on all channels. The 
action that accounts for the interchannel Coulomb interactions 
of same-chirality branches is 



i J = l '- 

(48) 

where (j>i, i = 1,2,3, are bosonic fields associated with 
each CS on one side of the gates. dx<t>i = —^^Pi{x) 
is the charge density of channel i and quantization is as- 



sumed based on the commutation relations [pi{x), Pj{x')] = 
— {ii'/2TT)SijdxS{x — x'). We can assume in general the intra- 
channel interactions Vu for i = 1, 2, 3 to be different, and the 
interchannel ones to fulfill Vij = Vji 7^ Vu. The diagonalized 
action can be written as 



S ■ 



/ da; dr 

1=1 



(49) 



where ipi = Mijipj and are eigenvalues of V. Since all 
channels are of the same chirality, the interaction-dependent 
AI matrix is unitary, {M'^M)ij — Sij, as opposed to the 
case corresponding to the presence of opposite-chirality chan- 
nels in which case M^JM = J. The matrix elements of 



J are (J)i 



Sij sign(j) and J contains at least one el- 



ement of opposite sign compared to the rest, corresponding 
to the opposite-chirality channels present. This will result in 
(Af^^Af)y ^ 6ij, thus affecting scaling exponents of tunnel- 
ing processes. The correlation functions between electrons of 



branch i, Rf ~ e'^'^'/'^, have the form 



and the scaling dimensions are di ~ ^^{M'^ M)ij /v ^ 
iiTelevant for = 1/3 and unaffected by the interactions. 

Another point to consider is the question of stability of the 
channel structure under interchannel tunneling processes me- 
diated by impurities along the edge. For repulsions Vij — V 
( <Vii = C/), the action is diagonalized in terms of two neutral 
and one charged field. At low energies the neutral fields tend 
to antilock charge densities in the consecutive channels lead- 
ing to an oscillatory charge density profile across the edge {x 
direction), in analogy to semi-classical results obtained in Ref. 
32. Interchannel tunneling processes lead to the chiral sine- 
Gordon model cos{\/STTipn/v), n — 1,2, for neutral fields. 
(In fact this does not depend on Vij.) The cos operator is rel- 
evant only for 0^ < IGTr,^"* a condition not realized in these 
systems since /3'^ — Stt/i/'^ and v — 1/3. 

The multichannel edges that appear in the insulating state at 
high I Vg I and high V offer a realization of models suggested in 
recent papers by Kim et al.f* where finite-frequency current- 
noise experiments measuring cross-current correlations could 
observe the fractional statistics of QPs participating in the tun- 
neling processes. The phase of the relative wave function of 
two QPs picked up upon their exchange can be observed in the 
proposed noise experiments if their exchange is mediated by a 
third channel. The systems studied in this paper provide such 
a channel naturally. Care, nevertheless, must be taken to local- 
ize the tunneling processes to a small region at, for instance, 
near the constriction as depicted schematically in Fig. [13] by 
placing additional gates to increase the distance of the edges 
along the rest of their lengths. 

As I Vg I is lowered we expect the 1 /3 IS channel to eventu- 
ally close with the remaining 2/3 IS to extend under the gate. 
At large bias, the same argument as in high \Vg\ for the con- 
ductance applies, with G — ve^ /h ~ — y2i/-2 ^ _ 2^3 
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FIG. 13: Schematic illustration at high | and large bias of the tun- 
neling processes; cross-current correlations, that can in principle be 
measured in finite-frequency current-noise experiments. Such pro- 
cesses are influenced by the charge and statistics of the quasiparticles 
that tunnel and can be used to detect their statistics (Ref.:34{l. 



At low V, however, Coulomb interactions between opposite- 
chirality channels as discussed in Ref. are expected to be 
strong leading to perfect transmission through the constriction 
and exhibiting the symmetry ly = 1/3 — 2/3 between 
the low and high \ Vg \ observed in the experiments of Ref. 6. 

One point that needs further comment here is the neces- 
sity of comparing this multichannel edge description with the 
one of the Hall systems of 2/3 bulk filling which is known 
to contain counterpropagating modes^^^ when electron cor- 
relations are taken into account. Based also on recent density 
functional theory (DFT) calculations on the edge structure in 
QH systems of integer filling, around an antidot^^ such modes 
do not seem to arise at the edges of 2/3 ISs formed at the inte- 
ger edge. Therefore our description, we believe, holds even if 
electron correlations are taken into account. 



B. Low Gate Voltage 

The experimental findings of suppression of backscattering 
between the top and bottom edges in QH systems with gate- 
created constrictions at small values of | | , whether at bulk 
filling 1/ — 1/3 or at integer filling — 1, according to our 
model, is the consequence of strong Coulomb interedge re- 
pulsions across the junction region. Features that can be best 
exploited for a consistency check of this proposal (of interest 
also in the broader context of observation of the LL model) 
are the ones that depend on the LL parameter K. The ones 
that we discuss here are two: a if -dependent oscillatory be- 
havior of I-V,^'^ and the formation of a particle-hole order 
parameter in the junction region. 

The interedge Coulomb interaction in systems of concern 
here is much stronger in the junction region than on other por- 
tions of the edge. Therefore, in these systems the interactions 
define a piecewise Fermi velocity vp along the length of the 
non-chiral LL, formed by the counterpropagating edges on the 
opposite sides of the junction, vp = Vp/K in the junction re- 
gion and Vp in the other portions of the edge. Discontinuities 
in the value of vp, on the other hand, lead to physically ob- 
servable consequences in transport experiments. Namely, the 
requirement of continuity of the current at all points along the 
LL leads to the phenomena of partial Andreev reflections^'^ 



of edge waves at the interaction inhomogeneity boundaries 
(x = ztw OT X = L/2 ± w, where L is the perimeter of 
the QH edge and w the length of each gate. Fig. [T4b . On 
the other hand, the finite length of the QH line junction in- 
troduces a time (or frequency o^lj ~ vf/2w) scale; the time 
for the edge waves reflected from the impurity to move to the 
junction boundary, where they incur further reflection, and 
come back again. The superposition at the point contact of 
these multiply reflected waves and the incoming ones can be 
tuned to constructive or destructive interference by, either a 
variation of strength of the Coulomb interactions (thus vp and 
cjlj) that can be tuned by changing the separation distance be- 
tween the left and right edges, or by a variation of the length 
of the gates, leading to oscillatory behavior of the frequency 
Lu ^ cjlj of the current transmitted through the constriction. 
The intensity and amplitude of these oscillations, depending 
on the amplitude of the multiple reflections of the edge waves, 
varies as a function of K, too. Oscillatory behavior in I-V 
was observed in fact in the experiments^ and'" although a de- 
tailed analysis of its frequency dependence on junction length 
A = 2w, interaction strength, and most importantly a distinc- 
tion from oscillatory behavior originating from other sources, 
unrelated with the presence of electron correlations in the sys- 
tems (like Fabry-Perot ones"*", etc.), for purposes of interest 
here remains to be done. 

The other point of physical importance that arises due to 
electronic correlations in the junction region, is the formation 
of the particle-hole order parameter there: 

0(x,i) i?°'-^(x,t)i?°'(-x,t) . (51) 

The origin of these phenomena is the drag effect induced by a 
propagating charge density on one side of the barrier to the 
other side to which it is coupled by Coulomb interactions. 
Formally this appears in the splitting of the bosonic fields into 
left-right movers that physically extend on both edges. The 
geometry assumed here is that of a single edge enclosing a 
2DEG in a QH plateau of filling factor Fig. [14] (bottom). 
The Hamiltonian as usual is described in terms of charge 
density fluctuations Ti{x^x') — p(x)U{x,x')p{x'), where 
X, x' are coordinates along the perimeter of the edge and 
U{x,x') = S^E/Sp{x)6p{x')\^, is the interaction kernel. In- 
teredge interactions at low |V^| are of Coulomb type, Vd{x) = 
e^/ {e[x^+al'^]^^'^) and are concentrated only on the region of 
the gates, x e {L — w,w) and x £ {L/2 — w,L/2 + w), where 
L is the edge perimeter length, a is a dimensionless constant 
of order one and serves as a short-distance cutoff parameter 
Interactions on the same side of the barrier, Vs{x), are also 
of Coulomb type. Here of importance is the exponent of the 

expected algebraic decay with time of (^O"^ {x,T)d{x,0)^, 

which we calculate numerically. To this end the Hamiltonian 
is diagonalized in momentum space with matrices M'^, M, 
such that M^JAI = J, where matrix J has elements differ- 
ent from zero only along its main diagonal, unity in absolute 
value but half of them being negative. The size of the matrices 
chosen in our numerics is of the order of hundreds. Using the 
identity of Ref.^ for the scaling dimensions d of the particle- 
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FIG. 14: Scaling dimensions of the particle-hole order-parameter 
operator of an inhomogeneous Luttinger liquid. The length of the 
system is L = 2400 / b ■ The fluctuations are relevant only in and 
just outside the regions of repulsive interactions. Interactions take 
place on both ends of the edge loop in a length of w = 200 Ib each. 



hole order parameter, we obtain 



1 47r r 1 



E 



i8sin(9fea;) sin((j,„x) 



M{qk,qi) [N{qi;0) - N{qi-T - r')] Aft(g,, g„,)}, (52) 

where qk = ink/ L with k integers k — 0, ±1, ±2 . . . (peri- 
odic boundary conditions were assumed) and N{q\T — t') = 
exp{ — Ag I r — r' I } /2 and Xq are energy eigenvalues (eigenval- 
ues of JH). Numerical results for the scaling dimensions of 
the particle-hole order parameter are shown in Fig. [14] (top). 

In cases of closed constrictions, calculations can be done 
similarly. It is expected^ that left-right electron tunneling 
in the closed constriction limit is dual to the top-bottom in- 
teredge tunneling of the open constriction limit. In this case 
correlations of constituent operators have to be evaluated on 
each separate side of the gate. In case the interactions are over 
an extended region, the X±{q), where ± refers to the left and 
right Hall system on the left and right of the gates, each of 
filling 1^1 and 1^2, are given by 



A±(«) 



Vl + V2 



viV2Vl{q) 



1/2 



(53) 



If i/j = 1^2 then \±{q) = X{q), with X{q) = iy[V^^{q) - 
V^{q)]^^^\q\. In general, however, calculations for eigenval- 
ues have to be done numerically. 

The particle-hole order parameter is a relevant operator in 
the interaction region and it persists to be so even in a small 
interval outside this region, see Fig. [14] The relevance of this 
operator depends on the strength of the interactions. Exper- 
iments probing the distribution of charge of waves incident 
on the line junction like the ones discussed in Ref. 43, which 
can lead to the observation of the particle-hole order in and 
around the gate region, can possibly measure K as well as its 
Vg dependence. 



VI. CONCLUSIONS 

In the current paper we have explored the reasons for the 
occurrence of the insulating behavior at high \Vg\ at a con- 
striction created by electrostatic gates in a 2DEG in which 
the integer QHE has been established. First we examined the 
nature of interactions between edges in the constriction re- 
gion. More precisely, we answered the question of how image 
charges appearing on the gates influence the interedge inter- 
actions. We find that at low |V^| (small 2DEG-gate distances) 
the interedge interactions are well screened. At large values of 
I Vg I a sign reversal of interactions does not occur; interactions 
remain repulsive in nature at any voltage of 2DEG subsys- 
tems. Here the analytic electrostatic calculations have been 
done exactly with no approximations. We then studied the 
structure of the gate-induced edge in the presence of mag- 
netic fields. Smooth edges have an instability toward for- 
mation of compressible and incompressible stripes, around 
points of fractional (or integer at higher) fillings. It turns 
out that the experimental geometry and values of gate volt- 
ages used in Ref. 6 lead to formation of wide ISs whose pres- 
ence as shown in Fig. [T2I can lead to enhanced top-bottom 
fractional-charge backscattering or, equivalently, to a suppres- 
sion of left-right electron tunneling through the constriction. 
Moreover, the G-V characteristics of constriction tunneling 
are expected to be similar to the one of bulk fillings 1/3, as 
observed in the experiments.*' The consistency of models pro- 
posed here can be verified, at high |V"g|, in addition to noise 
experiments commented on in the last section, also by equi- 
libration experiments similar to those of Alphenaar et a/.— 
One can inject current on one channel of the integer bulk edge 
by the use of additional metallic gates that separates only the 
outermost channel, for instance, bringing it in contact with a 
current source. One can measure its redistribution in time as 
well as its splitting as it passes through the point-contact re- 
gion. At low I V^|, properties of the emerging inhomogeneous 
LL (discontinuity in value of Fermi velocity) with its tunabil- 
ity of interaction strength by the voltage of the gate as well as 
the finite size of the junction region can be exploited in consis- 
tency verification experiments. The /-C-dependent oscillatory 
behavior of I-V, and the particle-hole order parameters were 
proposed as signatures to be observed in the experiments. 
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APPENDIX A: POTENTIAL CONFIGURATION OF 
CHARGED DIELECTRIC STRIPES 

1. General form of charge distributions 

A 2D potential can be represented as the imaginary part of 
an analytic function F{() with C, = x + iz. For the deriva- 
tive dF/dC, = dz ImF + \dx ImF holds and therefore the 
components of the electric field are 



Er = — Im 



dF 



dC 



Ez = -Re 



dF 



dC 



(Al) 



The boundary conditions of Sect. are applied to Im F 
(/)^''- (x, z). In the depleted region there is 



dlmF 

dz 



Re 



dF 



dC 



t{x) 



a < \x\ <b 



(A2) 



Since the function cj)^*^'-' ■ {x , z) is constant along the x axis out- 
side the above intervals (at z = 0), Ex{x, 0) — 0, i.e.. 



dlmF 

dx 



Im 



dF 
dC 



, < a or |x| > 6 . (A3) 



Now another function, -D(C), is introduced whose imag- 
inary part is known everywhere along the x axis. From 
(IA2l l. (IA3b one knows Rc[di^/dC] in part of the real x axis, 
a < |x| < b, and also Im[dF/dC] in the remaining part of 
the X axis. Therefore, if one multiplies this function with an- 
other, known one, which is purely imaginary in the intervals 
a < \x\ < b and real otherwise, then a new function f{() is 
built whose imaginary part one knows everywhere along the 
X axis. 

Such an auxiliary function is i[{x^ — a^)(6^ — x^)]^/^. 
Therefore, the new function whose imaginary part is known 
everywhere along the x axis, is given by 

m-m)^ , Dic) = [{e-a'){b'-e)]'/' . 

(A4) 

We can generate the function in the whole complex plane by 
using the Schwarz equality 



/(C) - - 

TT 



da; 



D{x)[dF/d(:]{x) 
x-C 



(A5) 





FIG. 15: Top: a plot of the potential (jf^'^'{x, z — 0) created by a 
dielectric ^ for certain values of a = 11b, b = WIb- Bottom: 
a plot of the potential (x, z = 0) = ImF{Q, where F{() is 
given in dAllb . The selected values of b, a, a, b are b — —81b, 
d — —31b, a = llB,b — 101b- 



In our case dF/d( obeys Eqs. ( IA2l i. (IA3l l. 

The integral iA5i just reproduces the numerator of its inte- 
grand on the real axis including a factor i. For some analytic 
t{x) the result can be expressed as 



fiC)=iT{C)D{0+T{C,a,b) 



(A6) 



where a, b) is also analytic and is to be found from the 
boundary conditions. These are (a) lm[T{x, a, b)] = every- 
where on the real axis and (b) f{C,) — > as |C| — > oo, which 
is imposed by (IA5b . For Tn{x) = Tnx" one obtains 



L ^ J 

.F„(C,a,5) =r„C"+^ ^ C;^{^) {ab/ey , (A7) 



where C^, ^ (t) is a Gegenbauer polynomial, cf. Eq. (8.930) of 
Ref. [l9i and [•] represents the integer part. For constant t{x) 

the function /(C) results in {C^^ {t) = 1, CP (t) = -t) 
f{C)^ro[iD{0 + e-h,{a' + b')] . 
The electric field is already gained from iA4i via iAli . 



dF_ 

dC 



r{C) 



(A8) 



which has the expected singular square-root behavior at every 
edge. These are the only singularities caused by D{Q, and 
since it appears in the denominator, it should also not have 
any additional zero in the complex plane. Together with the 
requirements imposed for its construction, the function -D(C) 
is (almost) fixed. 
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The function F introduced is obtained by the integration 
of (|A8t along a curve within the considered half-plane from 
a to These integration boundaries fix the boundary poten- 
tials outside the dielectric stripes to zero in the gate region and 
(/)'^''' {x, z = 0) = A{a, b) in the 2DEG regions. The potential 
constant A{a, b) introduced in Eqs. (|4|l and ^ comes out to 
be: 



A{a,b) = -Tb 



Eik')-Ul + ^]Kik') 



, (A9) 



where k = a/b and fc' = (1 - a^/b^)^^^. 

For the general case of Sect. |III]the potential created by 
the dielectric stripes, after integration of f{0/D{C),^^'^* is 
obtained to be 



T Im 



(6- 



a + by 



;Fiip,k) (AlO) 



.4[(5-a)(a-&)]i/2 
i[(6-a)(a-6)]i/2 (^E{^,k) ~ F{ip,~k) 



(b-C)(C-a)(C-&) 
C-a 



1/2 



(C-«)} 



Here k is the modulus of the elliptic functions, k = 

— —1/2 

[(6 — a)(a — b)/{b ~ a){a ~ fe)] , and (p is given hy ip — 

arcsin [{b — a){( — a)/{b — a){( — a)]^^^. [The singularities 
from the square-root term and from ip in (lAl 11 1 at C ^ a just 
cancel (Fig. [TsT l.l It is worthwhile to notice that the modu- 



lus k relates to k — Sf. / 

fc= (l-fc)/(l + fc). 



2. Polynomial charge distributions 

The charge density for ISs is given by the density profile 
n{x) ( fT3T l. see Sect.|lV] Approximating To{n{x) — n/), tq = 
27re/e, up to terms containing the A^th derivative of n{x) at x/ 
we arrive at different polynomial approximations for intervals 
a' < < b', where we denote by t/^'' the coefficients of a;' 
for a given order N. 

For = 1 the (second-order) approximation is given by 
(|40] i with Tq^-* —^ XfTQn'f and T2^^ = iron^/xy. The 
function F from (lAll) is obtained by the integration of ( IA8b 
with terms from ( IA7I ) and D{x) replaced by D'{x) — 
[{x^ - a'2)(6'^ - x^)Y/'^, resulting in 



T^^h' [C-k-i[{E' -E{xP,k')) (All) 

+ i(i + P)(£;'-^(tA,fc')) 

+ i(3fc'4 + 8P)(Js:'-^^(V,fc'))]} 



Here C, — C/b', k ~ a'/b' and ip,k' are given by the same 
equations as for (O, with a and b substituted by a' and b'. This 



leads to the potential difference ^'^^ (|4TI ). The approximation 
for N = 2 (fourth-order) is given in (|46] |. One can partially 
reuse the results for N = 1 using different coefficients . For 
F this leads to 



+ fe"{(C' - ~k') + i[c(C' + ^(1 + P))D'{b'Ob" 

+ 1^(1 + fc') (5^'^ + 8P) (-ft:' - F{^, fc'))] } , 

and the value for A^^^^ (a', b') reads 



+ 40^4 



I{b) /I [a) of Sect. HITBl bv the difference. 



Also, as for the second order in Sect. lIV fl difference and sum 
of the equations ( l47b yield 

= 2^/ + frf {2k'^K' - (1 + fc'^)^;') 
+ iri^'fo'S (4fc'2(l + - (3 + 2fc'2 + 3k'')E') , 

where in the sum was inserted and To^-* was taken from 



APPENDIX B: POTENTIAL CONFIGURATION OF 
METALLIC COMPONENTS 

1. Closed form for the general potential 

The general solution for (?!i°''(x2, 2/2) for the capacitor rep- 
resented in Fig. [3] (right), of AK periodicity along the X2 di- 
rection, reads <j)''^- ^ V2 + {Vi - ¥2)0 /2K + Voy2/K' + 
{Vi - V2)/{2K)U, where 

J'P^^Z_^ „;„u/„. T^/\ cos(a„X2) (Bl) 



n=l 



Oin svahia.nK-') 



with OLn — mi/2K and, separated in — f^^^ + f'^\ is 
given in ( IB2I ) and ( IB4l i. 

The first term containing products of trigonometric func- 
tions and an exponential can be written as 



ia„(C2+C) _ „ia„(C2-C) 




= _c.H^R.|n„!i^(if±§ ,(B3, 

with C2 = ^2 + it/2, and where, in general, depending on 
the arguments inside the logarithm, a constant 41f An may be 
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added. The second term containing products of trigonometric 
functions and a hyperbolic sine can be written as 



/, 



(2) 



(B4) 



TT 71 1 



°° 2 g-2Q„JC' 



-2a„K 



7 sin(a„C2) sin(a„C) 



This result for the symmetric case can be obtained in a sim- 
pler way, where the limit C — > If in (IB lb just leaves the odd n 
with an alternating sign. There the summation, after a few ma- 
nipulations involving an intermediate derivation by (^2, leads 
directly to elliptic functions of modulus k. 



Since here just a product of two sine functions occurs, the 
following formula^** can be used: 

E°° g^" sin 2nQ! sin 2n/3 i?i(a + /3, g) , sin(a + /3) 
7Z — In ^r-, — 7— In - 



n=l 



n(l - 



z?i(a-/3,g) 



i(a-/3)' 



identifying a = 71(2 /4:K and (3 — ttC/AK. Because k and 
therefore fc', JC, fC', and q — exp{—TTK'/K} are fixed, the 

variable q related to enters as q — > ^/q into the sum re- 
sulting in 



2K 



I In 



(B5) 



For — /^^' + Z^^-* representing ( IB II ) the logarithmic sine 
terms cancel and (/!'°' (C2, (2) obtains the form ( l22l l. 

L;m;f to f/ze Symmetric Case — The limit in (|22] | is per- 
formed by setting by C ^ -fC, simplifying the first argu- 
ments of the i?i functions. Subsequently various formulae of 
Sect. 8. 1 of Ref. 19, are used. The quantity ^ corresponds to 
the modulus fctr = 2^/k/ (1 + k). After argument translation 
the functions are rewritten in terms of elliptic functions of 
modulus kti-, which then is transformed to k. The subsequent 
argument doubhng and translation yields 



The dynamical limit for C K, i.e., h —b, a —a, 
results from Vi = ¥2- The dependence of the edge position b 
on (Vi — Vg)/(Ta) is shown in Fig.[T6l 



b/a 




0.5 1 1.5 
(Vi-Vg)/(xa) 



FIG. 16: Plot of edge position b (see Fig.[2} as a function of voltage 
Vi — 14 is shown. Lengths are measured in units of half of the gate 
width, a, which here is taken to be a = 50 Ib- 



2. Equations Used for Ey 




(B6) 



The expression ( |26] | for i?"^ (C2) is rephrased u sing the Ja- 



cobi zeta function Z(u), Eqs. (16.34.1,2) of Ref. |28|, and its 
relation Eq. (17.4.28) of Ref. 28 to elliptic integrals E and el- 
liptic functions. Their modulus is fctr = 2\/fc/ (1 + fc), which 
is transformed to fc for both function types leading to 



<(C2 



, - V2 



K' 



2tt 



Re 



(C2 + C)§ + 



cn 




dn(^,fc) 


sn 







- (C ^ -C) 



[Notice also that the arguments of the incomplete elliptic in- 
tegral of second kind E{-, •) depend on the convention: see 
Refs. [m andii] Then using Eq. (8.157) of Ref. 19 for half 



arguments and subsequently Eq. (8.151 2.) of Ref.[T^ making 
the selection of the sign depending on the sign of the elliptic 
functions noticing < C < 2K. With this and the argu- 
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ment shift for elliptic integrals, Eq. (17.47) of Ref . l28l for the 
special arguments ^2 = 0, 2K, one gets the expressions 



V1-V2 



+ k' 



1 -sn(C- Jf) 



n 
_, 1 



E 



2i5(am(lC))-C- 



1 + sri(C - K) 



E±(2K) 



K' 



V1-V2 



2k 



'dn{C-K)-k' 



TV 
_, 1 



dn(C -K)+k' 



2E (am(lC))-c|: 
l+8n{C-K 



1 - sn(C - K) 



(B8) 



,(B9) 



where the second argument of am, E, sn, dn is k. 

Starting from the expression for si (Ref. one obtains 
the shift P of Eg. ( |20] |. The modulus k of the transformation 
is defined in ( |21|) as I{b)/I{a); its explicit form is given in 
Sect. 1111 Al as k — /S'^. From that the right-hand side of 
the equation sn(C — K, k) ~ P/I{b) can be written as 



_P_ ^ [{a-b)ia-b)]'/^-[{a-b){a - b)] '''' 
-^(^) [{d-b){a-b)]'^^^ +[{d-b){a 



(BIO) 



which allows to express the following terms: 



k' 



1 T sii(C - K) 
1 ± sii(C - K) 



2[ia^-b^)ia±~b^)Y- 



Si 



(BID 



dn(C -K)-k' 



_Atl{C -K) + k' 
b-b - [{a - b)ia - b)]^ 



[{a-b){a-b)y 



Si 



■ (B12) 



In jBl 11 1 we have used the short notation a±, respectively for 
a, a, and b±, respectively for b, b. 
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